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摘 ， 要 : 针对 非 线性 混合 整数 规划 问题 ， 本 文采 用 非 固定 多 段 映射 罚 函 数 法 处 理 约束 条 件 、 用 混合 整数 编 
码 技术 处 理 连 续 变量 和 整数 变量 ， 并 在 基本 差分 进化 算法 中 加 入 一 种 新 型 的 是 组 合 变异 算 子 和 一 
种 指数 递增 交叉 算 子 ， 由 此 构造 出 了 一 种 求解 非 线性 混合 整数 规划 问题 的 改进 差分 进化 算法 。 实 
验 表 明 ， 所 提出 的 算法 全 局 收敛 速度 快 ， 精 度 高 ， 重 棒 性 强 。 
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1 3 引 启 


混合 整数 非 线性 规划 问题 (MINP) 广泛 应 用 于 机 械 、 化 工 、 资 源 管理 、 生 产 调度 、 生 物 、 
军事 等 领域 。 求 解 MINP 的 传统 方法 有 分 支 定 界 法 、 广 义 Benders 分 解法 等 。 由 于 同时 含有 实 
数 变 量 和 整数 变量 ，MINP 是 一 类 NP 完全 问题 ， 随 着 变量 维 数 的 增加 ， 计 算 量 会 急剧 增 大 ， 
从 而 使 这 些 方法 存在 很 大 的 局 限 。 目 前 进化 计算 方法 广泛 用 于 解决 约束 优化 问题 ， 且 被 证 明 很 
有 效 凡 4。 近年 来 ， 有 许多 学 者 将 进化 算法 ， 如 遗传 算法 (GA)H-3、 模 拟 退 火 算 法 (SA) 向 、 进 
化 规划 (EP) 铝 、 粒 子 群 优化 算法 (PSO)ts7 和 差分 进化 算法 (DE) 等 用 于 求解 MINP 问题 。 这 
些 方法 一 般 都 取得 了 满意 的 效果 器 。 差 分 进化 算法 是 Storn 和 Price 于 1995 年 共同 提出 的 一 种 
采用 浮 点 矢量 编码 ， 在 连续 空间 中 进行 启发 式 随机 搜索 的 优化 算法 回 。 差 分 进化 算法 作为 一 种 
性 能 卓越 的 优化 算法 正 受到 日 益 关 注 ， 其 应 用 领域 也 越 来 越 广 。 为 了 将 DE 用 于 MINP 问题 ， 
文献 对 D 孔 的 变异 操作 进行 了 改进 ， 提 出 了 一 种 对 于 整数 变量 直接 在 整数 空间 进行 优化 计 
算 的 改进 差分 进化 算法 。 文 献 [10] 对 变异 矢量 采用 四 含 五 入 方法 进行 取 整 ， 使 之 适合 于 0-1 整 
数 规划 和 整数 只 包括 0,1 的 优化 问题 。 对 于 其 它 MINP 问题 ， 同 时 进行 向 上 取 整 和 向 下 取 整 运 
算 ， 得 到 两 个 试验 向 量 ， 从 而 扩大 了 寻 优 空间 ， 有 利于 提高 算法 搜索 到 最 优 解 的 鲁 棒 性 ， 取 得 
了 满意 的 效果 。 

本 文 在 文献 [10] 的 基础 上 ， 采 用 非 固定 多 段 映 射 罚 函 数 法 处 理 问 题 的 约束 条 件 ， 用 混合 整 
数 编码 技术 处 理 连续 变量 和 整数 变量 ， 结 合 差 分 进化 算法 两 种 不 同 变异 方式 的 特点 ， 给 出 一 
种 线性 递减 加 权 因 子 的 凸 组 合 n-13 变异 方案 ， 并 采用 指数 递增 交叉 算 子 提高 算法 的 全 局 搜索 
能 力 和 收敛 速率 ， 用 四 个 经 典 测试 函数 对 改进 算法 进行 了 测试 ， 实 验 结果 表明 了 所 提出 的 改 
进 DE (Modified DE, MDE) 算法 用 于 求解 MINP 问题 的 有 效 性 。 
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2 ”MINP 问题 描述 
一 般 混合 整数 非 线性 规划 问题 可 描述 为 
min f(z)= (7x,Y) 
S.t. gi(z,Y) < 0， i 二 二 MALS 


(1) 
hj(z,y) =0,， j=1,2,..,), 


<r<r, yy<y<y, 


式 (1) 中 zx 为 连续 矢量 ,yy 为 整数 矢量 ，(z7,z0) 和 (yz,y0) 分 别 为 相应 决策 矢量 的 下 界 和 上 
界 ，gi(x,y) < 0 为 不 等 式 约束 条 件 ，h;(zx,y) = 0 为 等 式 约束 条 件 。 这 些 约束 条 件 通 常 都 是 非 
线性 的 ， 因 而 MINP 问题 一 般 难 以 求解 。 对 于 等 式 约束 ， 可 以 令 hj(z,y) < 0 和 hj(zx,y) > 0 化 
为 不 等 式 约束 来 处 理 。 通 常 所 构造 的 广义 目标 函数 具有 如 下 形式 


F(zx,Yy) 一 f(z, y) 起 6(t) H(z, Y), (2) 


其 中 了 (x) = (x,y) 表示 原 目标 函数 ，65(t)H(zx,y) 称 为 惩罚 项 ，6(t) 表示 惩罚 力度 ， 五 (z,y) 为 惩 
罚 因 子 。 式 (2) 中 如 果 在 求解 约束 优化 问题 (COP) 的 整个 过 程 中 6(t) 固定 不 变 ， 则 称 为 固定 罚 
函数 法 ， 相 反 则 称 为 非 固定 罚 函 数 法 。 国 定 罚 函数 法 也 称 为 精确 罚 函 数 法 ， 在 整个 求解 过 程 中 
只 需 计 算 一 个 目标 函数 ， 但 惩罚 力度 6(t) 难 以 选择 。 如 果 太 小 ， 则 收敛 减 慢 ， 甚 至 出 现 不 收敛 
的 情形 ， 如 果 太 大 ， 则 会 在 可 行 域 边界 使 目标 函数 产生 病态 ， 给 计算 带 来 困难 名。 而 非 固定 罚 
函数 法 则 在 求解 的 过 程 中 ，6(t) 随和 迭代 次 数 上 是 变化 的 ， 将 约束 问题 转化 为 一 系列 无 约束 问题 
来 求解 。Parsopoulos 和 Vrahatis 用 PSO 算法 求解 COP 时， 提出 了 一 种 称 之 为 非 固定 多 段 映 
射 罚 函数 法 219， 对 约束 条 件 进 行 处 理 。 本 文 在 文献 [10] 的 基础 上 ， 对 分 段 映射 函数 进行 了 
适当 修正 。 设 


m+2l 

H(z,9y) = > 0(gi(z,y)) gi(z, gr (3) 
$=1 

qi(z,y) = max{0, pi(7,Y)), (4) 


gi(T,Y), ? 一 1 2 ) 772 
Di = ji-m 人 (ZI 人) 一作 十 二 人 十 2 十 了 
hi_m_tz,Yy), i=m+tltlm+lit+2,. 7 十 2 


(qi(z,Y)) = | 1 者 gi(z,y) <1 
2， 否则 ， 
4， 著 qi(z,y) < 0.001, 
8， 车 0.001 < gi(z,y) < 0.1， 
12， 车 0.1<gi(zx,y) 1， 
24， 和 否则， 
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其 中 (4) 式 表 示 对 约束 的 违背 程度 ，(5) 式 表示 惩罚 函数 的 强度 ，(6) 式 为 分 段 映射 函数 ，(7) 式 
为 随和 欠 代 次 数 变化 的 惩罚 力度 。 这 样 根 据 约 束 违背 程度 的 大 小 ， 自 适应 选择 不 同 的 惩罚 项 ， 将 
约束 优化 问题 转化 为 一 系列 无 约束 问题 来 处 理 。 


3 ”改进 的 差分 进化 算法 


差分 进化 算法 采用 浮 点 数 编码 ， 在 连续 空间 进行 优化 计算 ， 是 一 种 求解 实数 变量 优化 问题 
的 有 效 方法 。 要 将 DE 用 于 求解 MINP 问题 ， 必 须 对 DE 进行 改进 。DE 的 基本 操作 包括 变异 、 
交叉 和 选择 操作 ， 与 其 它 进化 算法 一 样 也 是 依据 适应 值 大 小 进行 操作 。 根 据 DE 算法 的 特点 ， 
只 要 对 变异 操作 进行 改进 就 可 以 将 DE 用 于 MINP 问题 。 本 文中 对 变异 后 的 矢量 同时 向 上 和 向 
下 取 整 运算 ， 然 后 择优 选择 变异 个 体 进行 交叉 操作 ， 从 而 扩大 了 和 寻 优 空间 。 对 于 整数 变量 只 
取 0,1 两 个 变量 的 优化 何 题 ， 在 取 整 运算 时 采用 四 侈 五 入 法 。 

3.1 变量 描述 与 初始 化 

DE 由 N 个 矢量 (zt,yt) (i = 1,2,.… ,NN，t 表 示 第 t 代 ) 构 成 种 群 ， 这 NN 个 矢量 在 搜索 空 
间 进 行 并 行 直接 的 寻 优 。 设 实数 变量 的 维 数 为 n"， 整 数 变量 维 数 为 4， 则 (zt, 优 ) 可 表示 
为 (ZL 7Z2 ,znisyliyy2i,… ,ydi)。 初 始 化 时 ， 根 据 式 (8) 和 (9) 分 别 对 实数 变量 和 整数 
变量 进行 初始 化 。 对 于 整数 变量 ， 先 在 实数 空间 进行 随机 取 值 ， 然 后 同时 向 上 和 向 下 进行 取 整 
运算 ， 择 优选 择 其 中 一 个 个 体 进 入 初始 种 群 ，( 对 于 0-1 变量 ， 先 在 0,1 空间 进行 随机 取 值 ， 然 
后 四 舍 五 入 取 整 )， 得 到 对 应 的 整数 变量 


z= 72’ +rand( )* (7 ~ x’), (8) 
W= [y+rand( )* (y” —y")], (9) 


式 中 rand( ) 为 01 之 间 的 均匀 分 布 随机 数 ，[#] 表示 相应 的 取 整 运算 。 

3.2 ”变异 操作 

DE 最 基本 的 变异 成 分 是 父 代 的 差分 矢量 ， 每 个 矢量 对 包括 父 代 两 个 不 同 的 个 体 (zt1, Zra)。 
根据 变异 个 体 的 生成 方法 不 同 ， 形 成 了 多 种 不 同 的 差分 进化 算法 方案 门 ， 其 中 DE/rand/1 方 案 
全 局 搜索 能 力 强 ， 但 收敛 速度 慢 ，DE/best/1 方 案 局 部 搜索 能 力 强 ， 精 度 高 ， 收 敛 速度 快 ， 但 
会 加 大 算法 陷入 局 部 最 优点 的 可 能 性 。 结 合 这 两 种 不 同 变异 方式 的 特点 ， 本 文 用 一 种 新 的 变异 
方案 ， 其 变异 操作 方程 为 


(Zuo) 一 入 (Z73， ya) (d= 入 )(zbest， Ybest) +F* ((z71, yt1) Ey (zr2， yr2)), (10) 


式 中 (Zi1, Yr1)) (zr2, Yr2), (Zr3, Yr3) 为 互 不 相同 的 随机 个 体 ， (zbest， ybest) 为 种 群 中 适应 值 最 
好 的 个 体 ， e [0,2] 为 缩放 因子 ， 式 中 Xe [0,1]。 若 和 = 1， 则 (10) 式 等 价 DE/rand/1 方 案 ; 
若 和 = 0， 则 (10) 式 等 价 DE/best/1 方 案 。 一 般 来 说 ， 良 好 的 算法 要 求 在 初始 阶段 有 较 强 的 全 
局 搜索 能 力 ， 尽 可 能 发 现 多 的 可 能 全 局 最 优点 ， 而 在 后 阶段 则 应 有 较 强 的 局 部 搜索 能 力 ， 以 提 
高 算法 的 精度 和 收敛 速率 。 为 此 ， 将 DE/rand/1 方 案 和 DE/best/1 方 案 进行 凸 组 合 ， 引 入 线 
性 递减 加 权 因 子 和， 将 入 设置 为 退火 因子 ， 如 


入 = (Tinax — t)/Tinax (11) 


所 示 ， 其 中 Tmax 表示 最 大 迭代 次 数 ，t 表 示 当 前 迭代 次 数 。 和 在 搜索 过 程 中 由 1 逐渐 变化 为 0， 
使 得 zi 的 权重 逐渐 减 小 而 zs 的 权重 逐渐 增加 ， 从 而 保证 算法 既 有 较 强 的 全 局 搜索 又 有 较 快 
的 收敛 速率 和 搜索 精度 。 
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3.2.1 ”实数 变量 变异 操作 


Ty = Mt3 + (~ NThost + F(Tr1 ~ Zr2). (12) 
3.2.2 各 数 杰 量变 算 控 作 
对 于 整数 变量 ， 同 时 对 变异 矢量 向 上 和 向 下 取 整 ， 从 而 得 到 两 个 不 同 的 变异 矢量 。 然 后 对 
两 个 变异 个 体 择优 选择 ， 进 行 交 叉 操 作 ， 这 样 有 利于 扩大 搜索 范围 ， 从 而 有 利于 提高 算法 搜索 
到 全 局 最 优点 的 鲁 棒 性 。 对 于 0-1 整 数 规划 和 整数 只 为 0,1 的 情况 ， 则 用 四 舍 五 入 的 方法 进行 
取 整 。 对 整数 变量 进行 变异 操作 的 方程 为 


一 [Avts (1 Ws A) West 宇 人 (yii yz)], (13) 


其 中 [类] 表示 相应 的 取 整 运算 ， 通 过 对 变异 矢量 进行 改进 之 后 ，DE 就 可 以 用 于 实数 变量 ， 整 数 
变量 和 混合 整数 变量 优化 问题 。 

3.3 ”交叉 操作 

DE 利用 交叉 操作 以 保持 种 群 的 多 样 性 。 对 于 群体 中 第 i 个 个 体 (zi, yt)， 将 与 (7z,,%) 进 
行 交叉 操作 ， 产 生 试 验 个 体 (zz,yr)。 为 保证 个 体 (zt,yt) 的 进化 ， 首 先 通 过 随机 选择 ， 使 
得 (zr;,yzr) 至 少 有 一 位 由 (zu 如) 贡献 ， 而 对 于 其 它 位 ， 可 利用 一 个 交叉 概率 因子 CR， 决 
定 (zz7,yr) 中 哪 位 由 (zu, 加 ) 贡献 ， 哪 位 由 (zt, yt) 贡献 。 交 叉 操 作 的 方程 为 


Zzvj， 车 rand()<R 
zTj = j= 1,2,...,n, (14) 
zt;， 否则 
Yvj; 若 rand( ) < 及 
277 三 1 = 1,2 i ,d, (15) 
其 中 rand( ) 为 ,站 之 间 的 均匀 分 布 的 随机 数 ，CR & [0,1]。 由 
CR=CRn 二 (CR -CRain)eraG-tTnao) ，a=30， b=3, (16) 


可 知 ，CRR 越 大 ， 则 (zy 对 (zzr),yr) 的 贡献 越 多 ， 当 CR = 1 时 ，(zu 加 ) = (ZT,YyT)， 
有 利于 局 部 搜索 和 加 速 收敛 速率 ;CR 越 小 ， 则 (zc 的 ) 对 (zr,y7) 的 贡献 越 多 ， 当 CR = 
0 时 ，(zi yi) 对 (zz,yr) 有 利于 保持 种 群 的 多 样 性 和 全 局 搜索 。 由 此 可 见 ， 在 保持 种 群 多 样 性 
与 收敛 速率 之 间 是 矛盾 的 。 和 良好 的 搜索 策略 应 该 是 在 搜索 的 初始 阶段 保持 种 群 多 样 性 ， 2 
全 局 搜索 ， 尽 可 能 得 到 多 个 可 能 全 局 最 优 的 个 体 ， 而 在 搜索 的 后 期 应 加 强 局 部 搜索 能 力 ， 
提高 算法 的 精度 。 基 于 这 种 思想 ， 本 文采 用 指数 递增 交叉 概率 因子 CR 的 方法 ， ee 
次 数 的 增加 而 由 小 变 大 ， 初 始 阶段 (zt 倪 ) 对 (27,yr) 贡献 多 ， 提 高 全 局 搜索 能 力 ， 而 在 后 期 
则 (zi 开 ) 对 (zz,yr) 页 献 多 ， 提 高 局 部 搜索 能 力 。 设 CRmin 为 最 小 交 又 概率 ，C Rmax 为 最 大 
交叉 概率 ，t 为 当前 迭代 次 数 ，Tihmax 为 最 大 迭代 次 数 ， 则 CR 由 (16) 式 确定 。 

3.4 选择 操作 

DE 采用 “ 贪 栖 ”的 搜索 策略 ,经 过 变异 与 交叉 操作 后 生成 的 试验 个 体 (zT,yr) 与 (zt 罗 ) 进 
行 竞争 ， 只 有 当 (zzr; 条 ) 的 适应 度 值 较 (xt, 孝 ) 更 优 时 才 被 选 作 子 代 。 否 则 ， 直 接 将 (zf 2) 作 
为 子 代 。 选 择 操 作 的 方程 为 


(zt ttl yt tt1) = (ZT， yT), 车 F(zT, YT) < 下 (zi yt), 
(zi, 如)， 否则 . 


(17) 
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3.5 ”算法 描述 

综合 以 上 对 DE 的 改进 ， 我 们 提出 求解 MINP 的 改进 DE 算法 的 流程 如 下 : 

步骤 1 初始 化 种 群 规 模 N， 一 般 为 参数 维 数 的 10 倍 ， 收 缩 因子 已 ， 交 叉 概 率 上 下 
界 CRmax 和 CRmin， 控 制 误差 e， 在 每 个 变量 的 定义 域内 按 式 (8) 和 (9) 随 机 初始 化 每 一 个 
个 体 。 设 置 差分 进化 算法 最 大 迭代 次 数 Tnsx， 置 当前 迭代 计数 器 + = 0， 置 罚 函 数 当 前 迭代 计 
数 器 t = 0; 

步骤 2 按 (4),(5) 和 (6) 式 计算 每 个 个 体 每 个 约束 条 件 的 惩罚 因子 ; 

步骤 3 按 (3) 式 计算 每 个 个 体 的 所 有 约束 条 件 的 惩罚 因子 H(zx); 

步骤 4 按 (2) 式 计算 每 个 个 体 的 适应 值 ， 求 出 最 优 适 应 值 及 最 优 个 体 xt; 

步骤 5 判断 6(t)H(z,y) 是 否 达 到 精度 要 求 ， 若 是 ， 则 退出 。 否 则 差分 进化 算法 开始 从 
代 ， 执 行 下 一 步 ; 

步骤 6 在 种 群 中 随机 选择 三 个 与 (zi gt) = 1,2,… , 入) 不同 的 个 体 ， 按 (12), (13) 式 进 
行 变异 操作 ， 生 成 变异 个 体 (x,,y,); 

步骤 7 按 (14),(15) 和 (16) 式 进行 交叉 操作 ， 生 成 试验 个 体 (zx7, yr); 

步骤 8 按 (17) 式 进行 选择 操作 ， 生 成 t 十 1 代 个 体 (zit1, gt 

步骤 9 +=t+1， 返 回 步骤 2。 


4 ”数值 例子 


为 验证 本 文 算法 求解 整数 规划 和 混合 整数 非 线 性 规划 问题 的 有 效 性 ， 我 们 采用 如 下 四 个 典 
型 测试 函数 进行 实验 [215161， 
例 1 
min fi(x) = (zl 十 107z3)2 十 5(zs 一 Z4)2 十 (za ~ 223) + 10(7x1 一 10z4)4 
| st. zi€ [~100,100], i = 1,2,3,4, 
已 知 最 优 解 (zi, x2, za Za, 了 *) = (0,0,0,0,0)。 
例 2 
min fo(z,y)=27+Yy 
st. 1.25 一 Z2 一 22 <0， 
r+y<z1.6, 
0<Z<1.6， 
y € {0,1}, 
已 知 最 优 解 (x,y, f*) = (0.5, 1, 2)。 
例 3 
min fa(z,y) = 1 +(y 12+(y -1 -ln(y+1) 
二 (zl 一 1 十 (zl 一 2)2 十 (za 一 3)?2 
st 们 十 加 十 四 十 Z1 十 zz 十 T3 葡 5， 月 十 23 < 4.64， 
V3 十 ZT? 十 ZZ 十 3 <5.5, Yrz2 <1.8, ys+itrzs < 2.5, 
Vit+ri <1.2, y+r2 <1.64, yet+z1 <1.2, RR+7x3 < 4.25, 
Ti1, T2, 23 > 0, Yi, yo2, y3, ys € {0,1}, 
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已 知 最 优 解 (zt za Za, ,2,Y8;y4; 产 ) = (0.2, 1.28061, 1.95448, 1,0,0,1,3.557463)。 
例 4 
min fa(z,y) = —5.357854(z1)? -0.835689y zs — 37.29329y, -+ 40792.141 
St. 41 十 a2y2T3 + QA3yiT2 一 Q47173 < 92, 
Q5 + a6yazZ3 + aryiy2 十 as22 — 90 < 20, 
Qg 十 Qi0T17T3 十 Q11V171 十 Q12Z122 一 20 < 5， 
27 < 21,22);73 < 45， 
yn € {78,... ,102}, va € {33,... ,45}, 


已 知 最 优 解 (zl zz, za 加, 产 ) = (27,35,8175,27,78,36,3.22174 x 104)， 至 于 各 个 参数 的 取 
值 见 表 1。 


表 1: 例 4 中 各 参数 的 取 值 
















0.0022053 
0.0021213 


85.334407 Q2 0.0056858 0.0006262 


80.51249 | os | 0.0071317 | ay | 0.0029955 | as 
9.300961 0.0047026 0.0012547 | a12 


实验 中 种 群 规模 NN 取 变 量 维 数 的 10 倍 ，F = 0.5， CRwin = 0.5，CRwax = 0.9， 例 1 至 
例 3 的 最 大 迭代 次 数 为 50， 例 4 的 最 大 迭代 次 数 150。 为 减 小 随机 干扰 ， 每 一 问题 都 重复 20 次 
实验 。 表 2 为 每 一 问题 的 最 优 解 ， 平 均 最 优 解 及 标准 差 的 统计 结果 。 


表 2: 20 次 实验 结果 统计 表 









































澡 数 已 知 最 优 解 平均 值 标准 差 
例 1 | 0 0 0 
例 2 2 2 | 0 
例 3 | 3.557463 2.6368528 2.6368528 | 0 
例 4 | 3.221742e 十 004 | 3.221742e 十 004 | 3.221742e+004 | 0 





从 表 2 的 最 优 解 和 平均 值 可 以 看 出 ， 改 进 的 DE 算法 在 求解 问题 fs 时 最 优 值 明显 好 于 
文献 [10,17]， 对 于 文献 [10] 的 最 优 值 为 3.557463， 而 本 文 提 到 的 改进 算法 所 求 的 最 优 值 
为 2.6368528。 从 标准 差 可 以 看 出 ， 改 进 的 DE 算法 全 局 收敛 的 鲁 棒 性 好 。 对 于 问题 fi ~ 方 ， 
只 需 经 过 50 次 迭代 运算 即 可 得 到 最 优 解 ， 而 对 于 问题 户 也 只 需 150 次 迁 代 次 数 ， 说 明 算 法 收 
敛 速度 快 。 图 1 为 改进 DE 算法 求解 上 述 四 个 问题 时 20 次 平均 最 优 解 的 进化 曲线 。 


5 ”结论 


采用 非 固 定 多 段 映射 罚 函 数 法 处 理 问 题 的 约束 ，DEE 算 法 的 变异 矢量 同时 向 上 和 向 下 进行 取 
整 运算 ， 而 对 整数 只 包括 {0,1} 的 混合 整数 非 线性 规划 问题 则 用 四 舍 五 入 方法 进行 取 整 运算 ， 
使 改进 的 DE 算法 适应 于 求解 整数 规划 和 混合 整数 非 线 性 规划 问题 。 为 保证 种 群 的 多 样 性 和 提 
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高 算法 的 收敛 速度 ， 结 合 差分 进化 算法 两 种 不 同 变异 方式 的 特点 ， 引 入 线性 递减 加 权 因 子 的 丁 
组 合 变 异 策略 ， 并 采用 指数 递增 交叉 算 子 ， 对 四 个 典型 整数 规划 和 混合 非 线性 整数 规划 问题 进 
行 了 测试 。 实 验 结果 表明 ， 改 进 的 DE 算法 收敛 速度 快 ， 精 度 高 ， 全 局 搜索 鲁 棒 性 好 ， 是 一 种 
求解 混合 整数 非 线性 规划 问题 的 有 效 方法 ， 可 应 用 于 解决 各 种 实际 工程 优化 问题 。 





X 108 
1 一 2.14 厂 





Average best finess value 
站 
Average best fitness value 




















' Cn 
38 | ET 2 
f 


0 
3 


Avarage best finess value 
Average best fitness value 


© 























3 
| 
gL 
b 
所 


图 1: 改进 DB 算法 求解 例 1 至 例 4 时 20 次 平均 最 优 解 的 进化 曲线 
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Improved Differential Evolution Algorithm for Nonlinear 
Mixed-integer Programming Problems 
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Abstract: For nonlinear mixed-integer programming problems, we use the non-stationary multi-stage 
mapping penalty function to handle constraints and use mixed integer coding technology to handle real- 
valued variables and integer-valued variables. At the same time, a new convex combination mutation 
operator and an exponent increased crossover probability operator are integrated into the differential 
evolution algorithm. We construct an improved differential evolution algorithm for nonlinear mixed- 
integer programming problems. The experiments show that the proposed algorithm has fast global 
convergence speed, high computational precision and good robustness. 

Keywords: global optimization; nonlinear mixed-integer programming; non-stationary multi-stage 


mapping penalty function; differential evolution algorithm 
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